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We calculate the bound state properties of stoponium using lattice formulation of nonrelativistic 
effective field theory for stop which is moving nonrelativistically in the rest frame of stoponium. Our 
calculation method is similar to that employed in lattice nonrelativistic quantum chromodynamics 
(NRQCD) studies for charmonium and bottomonium. Using 16® x 256 quenched lattice gauge 
held conhgurations at a~^ = 50(l)GeV, we obtain the stoponium mass and the lattice matrix 
element which is related to the wavefunction at the origin for the IS' state and Hnd that the lattice 
|Ris(0)|®/Mj_g is 3.5 ~ 4 larger than that from a potential model calculation for 200GeV < Mis < 
SOOGeV. 


I. INTRODUCTION 

After the discovery of Higgs particle by ATLAS [T] 
and CMS at LHC, detailed measurements of its prop¬ 
erty become very important and urgent. These preci¬ 
sion measurements will open up a new opportunity to 
search for physics beyond Standard Model. In particu¬ 
lar, heavy particles which can decay into Higgs boson are 
under active experimental and theoretical investigations. 
Scalar top quark (stop), the supersymmetric partner of 
top quark and the Next-to-Lightest Supersymmetric Par¬ 
ticle (NLSP), is one of such possibilities. Then, stoponia, 
bound states (binding through SU (3) color gauge interac¬ 
tion) of stop and anti-stop may become interesting pro¬ 
vided that stop is long-lived enough to form a bound 
state. They can serve as a probe to stop, and decay of 
stoponium via the stop-anti stop pair annihilation into 
diboson states may be observed due to its distinct signa¬ 
ture PHZ]. 

With regard to the production cross section calcula¬ 
tion of stoponium, a next-to-leading order computation 
of perturbative part of the production cross section is 
available BE] and a resummed next-to-next-to-leading 
logarithm calculation is performed [10]. Turning these 
perturbative calculations into phenomenological compar¬ 
isons requires matrix elements for stoponium. With the 
observation that heavy quark moves slowly in the rest 
frame of quarkonium (u^ ^ 0.1 in bottomonium and 
u® ~ 0.3 in charmonium where v is the heavy quark ve¬ 
locity in the rest frame of quarkonium), nonrelativistic 
quantum chromodynamics (NRQCD) is developed and 
the quarkonium production cross section is given in terms 
of perturbative parts and nonperturbative matrix ele¬ 
ments mi. Similarly, the parton-level differential cross 
section for the stoponium (^') production in a collider 
can be given as 

da{ab -> 'k -|- A) = ^ da{ab —>■ ii[n] X){0'^[n]), (1) 

n 

where a, b are partons, t, t are stop and anti stop, n de¬ 
notes the angular momentum of the stoponium states 
and are generic forms of nonrelativistic stoponium 

production operators |5] by considering nonrelativistic ef¬ 


fective field theory (NREFT) for stoponium system since 
in the rest frame of stoponium for M ~ 0(100) GeV, 
is expected to be ~ 0(0.01) and as{Mv) ^ 0(0.1). 

Following m, one can relate the production matrix 
elements to decay matrix elements by crossing relation in 
leading order of . The decay matrix elements then can 
be related to nonrelativistic wavefunctions in the vacuum 
saturation approximation upto 0(u'^). For the 15 state 
stoponium, 

(O'") (0|xV|15)(15|V^t^|o) = |(0|xtV;|15)p, (2) 

when the leading u^-order term in the factorization ex¬ 
pansion is considered and a stoponium state is assumed 
to dominate in the intermediate states mi. Here, if de¬ 
notes nonrelativistic stop, x^ denotes nonrelativistic anti 
stop. This matrix element is related to the wavefunction 
at the origin m as 

A 

1(01x^115)^:..—|i?i5(0)P (3) 

and has been calculated on lattice for charmonium and 
bottomonium decay^ |12lH4] (see (15] for an improved 
lattice calculation for bottomonium system). 

So far, potential model estimates (e.g., [TB]) have been 
used for the stoponium masses and the stoponium wave- 
functions where typical bound states properties are sum¬ 
marized in with the mass of the S-wave stopo¬ 

nium state, Ms, and |i?5(0)| the value of stoponium ra¬ 
dial wavefunction at the origin. In , uncertainty 

from the wavefunction at the origin, |i?5(0)|, is more 
dominant than that from the mass. The mass of sto¬ 
ponium will be mostly from twice of the “free” stop mass 
and the “binding energy” will be just a few percent in 
stoponium mass and the uncertainty in binding energy 
will lead to sub-percent level uncertainty in stoponium 
mass. On the other hand, |i?5(0)| depends very much on 
the functional form of potential models naiiaiiHi. 


^ Note that in Eq. ^ the factor ^ instead of 3^; 
spin-less nature of stop 
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Clearly, using potential model estimates for the bind¬ 
ing interaction of stoponium is unsatisfactory since it 
introduces model-dependence which can not be system¬ 
atically improved, and makes the stoponium cross sec¬ 
tion calculation unsystematic as a whole, despite the 
improved perturbative calculations. Furthermore, po¬ 
tential models have a difficulty in obtaining right decay 
widths for a given state although they do better for rel¬ 
ative ratios for different states mmsi. Employing dif¬ 
ferent functional forms for potential models can lead to 
a large differences. For example, a potential model esti¬ 
mate for I i?s (0)1 defies a naive expectation that large stop 
mass (M > 100 GeV) would result in Coulombic behav¬ 
ior of the wavefunction, and exhibits substantial depar¬ 
ture from the Coulombic value of |i?5(0)| even at M ~ 1 
TeV [IH]. This suggests that there can be sizeable non- 
perturbative contribution to the bound state properties. 
Indeed, next-to-next-to-next-to-leading order calculation 
together with a scheme choice which is less sensitive to 
the long distance effect of QCD is necessary to under¬ 
stand the threshold behavior of the top anti-top S-wave 
pair production cross section m- 

In this work, we use lattice formulation of NREFT 
for stoponium and calculate the stoponium mass and a 
stoponium matrix element, |(O|x^'0|'I')p. Unlike a po¬ 
tential model calculation, lattice NREFT is based on the 
first principles of quantum field theory and allows sys¬ 
tematic study of errors associated with a given result. 
This effective lattice theory is similar to the lattice ver¬ 
sion of NRQCD [T^ which allows highly successful 
understanding of nonperturbative quarkonium physics 
(see e.g., [H]) except for the fact that stop is a spin¬ 
less particle. Our calculation for the stoponium property 
is performed on x N-r = 16^ x 256 lattices generated 
with “quenched approximation” and the lattice spacing, 
a~^ ~ 50 GeV for the stop mass range 1 < Ma < 8 (i.e, 
50 < M < 400 GeV). A large N^. is necessary to avoid 
deconfining effects. In the rest frame of stoponium, stop 
is expected to move slowly with velocity u(^ 1) and the 
size of stoponium should be smaller than those of typical 




As for the Monte Garlo data of SU{3) lattice gauge 
fields which are used in the calculation, they are gener¬ 
ated on 16^ X 256 lattice in “quenched approximation” 
at lattice bare coupling /3 = ^ = 8.751 using a single 
plaquette Wilson action. Multi-hit Metropolis algorithm 
together with interleaving over-relaxation algorithm m 
is used for the gauge field update. Each configuration 
is separated by 1000 Monte Garlo sweeps. To convert 
a lattice result into a quantity in a physical unit, one 
needs a lattice spacing as a function of the bare coupling 


quarkonia (rough estimate for the size of stoponium may 
be given by the self-consistency relation, v ~ as{Mv) 
with the size of r ~ {Mv)~^ and the stop mass M) and 
thus small Ng does not cause a significant finite volume 
effect. Since the momentum scale larger than the heavy 
particle mass is “integrated out” in NREFT and Ma is 
chosen to be ~ 1 , we need to consider only the lattice 
spacing scale in lattice NREFT for stoponium. 

We find that with 0{v^) NREFT Lagrangian, the lat¬ 
tice result for \Ris{0)\‘^/Mfg is factor 3.5 ~ 4 larger than 
a potential model estimate in [TB] for 200GeV < Mig < 
800GeV. Although further study is necessary (in partic¬ 
ular in view of the difficulty associated with quantifying 
the systematic uncertainty from the quenched approxi¬ 
mation), this implies that the stoponium production rate 
at LHG may be larger than the current estimates based 
on potential model. In the following, we briefly summa¬ 
rize the lattice method used in the calculation of stopo¬ 
nium properties (section [IT]) . Then, we present our main 
result in Sec. m and conclude with Sec. CYl 

II. METHOD 

The effective Lagrangian for nonrelativistic stop in 
leading order of is given by 

where ijj is a complex scalar field for stop which trans¬ 
forms as a SU{3) vector and y is that for anti-stop and 
Dr (D) are gauge covariant temporal (spatial) derivative 
under the strong interaction, SU{3). Note that this lead¬ 
ing order NREFT Lagrangian does not differ from that 
for NRQGD. The difference between stoponium NREFT 
Lagrangian and NRQGD Lagrangian occurs only when 
next-to-leading order NREFT Lagrangian is considered 
since the spin interaction (e.g., cr-B) in NRQGD is 
Additional terms for NREFT Lagrangian in can be 
systematically studied by including 


(D • E - E • D) V- + (D • E - E • D) x] . 


constant. In the scaling limit {Nf = 0 due to quenched 
approximation), 

aA^° = exp(- ^)ibo9T ^ = fig) ^ ^ 

(5) 

where bo = and ^(j^)^- We use the light 

hadron spectrum calculation in [51] to set the lattice scale 
since the experimental value for IP — l^" level splitting 
used in usual quenched lattice quarkonium calculations 
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is not available for stoponium. In m, at 6/g^ = 6.5, 
mp{mq —)■ 0)a = 0.200(4). We obtain a~^ = 3.84(8) 
GeV (Ai° = 5.10(1) MeV) from the lattice rripa and 
TOp(physical) = 768.1(5) MeV. This scale setting in¬ 
troduces ~ 2% systematic erroi|^ Thus, the scale for 
P = 8.751 is set as a"! = 50(l)GeV = 0.0039(l)fm. 

Under the background of these lattice gauge field con¬ 
figurations, nonrelativistic stop correlators for 0{v^) are 
calculated using the evolution equation, 

G(x,to) =S'(x), 

G(x,Ti) =^1 - U](x,T,_i) ^1 - G(x,Ti_i), 

( 6 ) 

where S{'k) denotes an appropriate complex valued ran¬ 
dom point source, diagonal in SU{S) color (random 
source improves the signal to noise ratio), and Hq is 
the lattice Hamiltonian corresponding to Eq. Q and 
G 4 (x,ri_i) is the time directional gauge field. The pa¬ 
rameter, k is introduced to stabilize large momentum be¬ 
havior of the lattice discretized evolution equation (see 
Table |T] for the choice of k). The gauge link variables are 
divided by “tadpole factor”, uq, which is chosen to be 

uo = (0|^TrUpiaq|0)5 (7) 

where Upiaq is a plaquette. For 0(v^) Lagrangian, a mod¬ 
ified evolution equation which includes the Hamiltonian 
for Eq. together with the improvement term for finite 
lattice spacing [20] can be used such as 

G(x,ro) =5'(x), 

x(l-5H)G(x,r,_i), (i>2), (8) 

where SH denotes the lattice Hamiltonian mentioned 
in the above. 

The zero-momentum stoponium correlators are then 
formed by combining the Green function for stop and 
that for anti stop and are summed over the spatial lattice 
sites. From the stoponium correlators, matrix elements 
are obtained by fitting 

Gs{r) 

n 

= + ■ ■ ■ . (9) 


^ however, this small error bar should be taken with caution. The 
quenched approximation effect in the hadron spectrum can be 
as large as 20% and the ambiguity due to the quenched approx¬ 
imation in the scale setting is far larger, (see e.g. m for the 
question of the scale setting in bottonionium system) 


For the lattice determination of the wave-function at the 
origin, |i?(0)p, the matrix element obtained from fitting 
Eq. is related to the nonrelativistic Coulomb-gauge 
fixed wave-function 


Ho = |(0|xV|l,5)r (10) 

in the leading order of with the normalization conven¬ 
tion for the radial wavefunction, dr \R{r)p = 1 

In NREFT formulation, the energy scale in the spec¬ 
trum is not known. In order to determine the mass of a 
state from the energy of state for a given channel from 
Eq. “energy shift” needs to be determined, 

Mpf = 2{ZmM — Eq) + Epf (11) 


where the mass of state, (V denotes the quantum 
number of the state), is given in term of the energy of 
the state (£"„), “zero point energy (Eq)”, and the mass 
renormalization (ZmM) |23| (or one can use the kinetic 
mass of nonrelativistic dispersion relation to determine 
the mass of a state [22]). In lattice NRQCD study of 
charmonium and bottomonium, this energy shift is de¬ 
termined either by the well measured experimental mass 
for one of the states (such as J/tjj and T) or by lattice 
perturbation theory computation |24| . Stoponium is not 
discovered yet and the perturbative estimate for the en¬ 
ergy shift using the lattice action used in this work is not 
calculated. Thus, we use tadpole improved mean field 
estimate for the energy shift in leading order |^, where 


Eq = —a ^ In 


uo 1 - 


aho 

Ijt 


2k 


with 


Zm = u 


-1 

0 



( 12 ) 


n is the integer parameter which is introduced in Eq. ^ 
to prevent the high momentum instability of the evo¬ 
lution equation. At /3 = 8.751(g^ = 1.4585), Monte 
Carlo simulation gives uq = 0.930049(1). Perturbatively, 

uq ~ 1-0.083^2- |2D] and uo « 0.94309 at ^ = 8.751. 

The difference between the lattice value of uq and the 
perturbative value of uq is ~ 1.4%, which suggests order 
of magnitude for nonperturbative effect in tadpole fac¬ 
tor, Uq at P = 8.751. Table [T| summarizes the mean field 
estimates of the energy shift for various Ma and n used 
for stoponium correlator calculation. 

There are various sources for systematic errors with 
nonrelativistic lattice calculation of stoponium proper¬ 
ties. An analysis similar to that in lattice NRQCD com¬ 
putation of quarkonium |20| can be applied: (1) rela¬ 
tivistic correction, (2) finite lattice spacing, (3) radia¬ 
tive correction, (4) finite lattice volume effect, and (5) 
light-quark vacuum polarization. In this work, there are 
an additional sources of errors. First is the determina¬ 
tion of the zero energy shift: due to the rest mass of 






4 


Ma 

k 

aho 

Zm 

Eoa 

2[Zm^ci — Eoa) 

1.0 

4 

0.209385 

1.04700 

0.285172 

1.52367 

2.0 

2 

0.104926 

1.04700 

0.178845 

3.83034 

2.5 

1 

0.083941 

1.03008 

0.158271 

4.83339 

4.0 

1 

0.052463 

1.04700 

0.125682 

8.12470 

5.0 

1 

0.041971 

1.05265 

0.114935 

10.2966 

6.0 

1 

0.034975 

1.05641 

0.107803 

12.4613 

7.0 

1 

0.029979 

1.05910 

0.102724 

14.6219 

8.0 

1 

0.026232 

1.06111 

0.098923 

16.7799 


stop mass Ma, where the vertical axis is in the loga¬ 
rithmic scale. By fitting these stoponium correlators, we 
obtain the energy of IS* state, Eiga and the amplitude, 
Table i summarizes the fit results for each stop 
mass with the fit range 60 < r/a < 100 where the fit 
range was chosen by locating the plateau region of the 
effective mass plot (Fig. [^. The error bar is from sin¬ 
gle elimination Jackknife error analysis of the fitted Eisa 
and From these lattice quantities, Eisa and A^a?, 

we get Mis and |i?is(0)p by use of Eq. (Ill and Eq. 
(101 (the last two columns in Table 0. 


TABLE I. mean field estimates of the zero energy shift 


the particle and its renormalization, NREFT has a un¬ 
determined “zero energy shift”, which can be fixed by an 
experimental mass of one of quarkonium states or by a 
perturbative calculation. Our use of a mean field theory 
approximation for the shift introduces a source of 
systematic error. Additional error concerns a bare stop 
mass used in the lattice calculation. An accurate calcu¬ 
lation of the stoponium masses requires “tuning of stop 
mass” which compares a stoponium kinetic mass in non- 
relativistic dispersion relation for a given stop mass with 
an experimental stoponium mass and changes stop mass 
until the kinetic mass is equal to the experimental stopo¬ 
nium mass. Since stoponium is not discovered yet, this 
tuning procedure can not be performed and the lattice 
result in the following contains a uncertainty from the 
imprecise tuning of stop mass. 


III. RESULT 


16^256, P = 8.751 



FIG. 1. The lattice non-relativistic correlators for the S-wave 
channel with 0{v^) lagrangian on 16® x 256 lattices 

Fig. 0 shows S-wave stoponium correlators, G(r), 
on 16® X 256 lattice volume calculated with the u® La¬ 
grangian (Eq. ([^) and the evolution Eq. ([^ for each 


Ma 

k 

Eis a 

Ao a® 

Mis (GeV) 

mor/Mfs 

1.0 

4 

0.1619(2) 

0.00507(4) 

84.28(1) 

4.43(5) X 10"® 

2.0 

2 

0.1688(1) 

0.02377(6) 

200.0(1) 

1.557(7) X 10“® 

2.5 

1 

0.1671(1) 

0.03599(7) 

250.0(1) 

1.205(5) X 10"® 

4.0 

1 

0.1553(1) 

0.1237(2) 

414.0(1) 

9.13(3) X 10-^ 

5.0 

1 

0.1455(1) 

0.2353(3) 

522.1(1) 

8.66(3) X 10-^ 

6.0 

1 

0.1344(1) 

0.4014(4) 

629.8(1) 

8.41(3) X 10-^ 

7.0 

1 

0.1220(1) 

0.6173(5) 

737.2(1) 

8.07(3) X 10-^ 

8.0 

1 

0.1087(1) 

0.8613(6) 

844.4(1) 

7.49(3) X 10-^ 


TABLE 11. Eis and Ao from lattice calculation (in lattice 
unit) on 16® x 256 lattices with NREFT Lagrangian. The 
result is based on 400 stoponium correlators and the error bar 
is from the jackknife analysis of the 1-exponential ht 


Figure|^shows lattice |i?is(0)|®/M®g (the column 6 of 
Table III as a function of lattice Mis (the column 5 of 
Table III in the range of lOOGeV < Mis ^ SOOGeV. In 
the figure, a line for a potential model result from m is 
also drawn for a comparison, where the A = 300 MeV pa¬ 
rameterization for Mis and \Ris{0)\'^/Mfg is used. The 
figure shows that the result from lattice NREFT calcu¬ 
lation is larger by factor ~ 4 at Mis ^ 200 GeV and 
by ^ 3.5 at Mis ^ 800 GeV than that from a potential 
model calculation. 

Let us consider the magnitude of the systematic er¬ 
rors in our lattice calculation of stoponium. First we 
consider the finite spacetime volume effect in Table 0 
by comparing with the results from two different lat- 
tice vo lumes. Table III from 20® x 256 lattices and Ta¬ 
ble IV from 12® x 256 lattices (both at /3 = 8.751 with 
the v"‘ Lagrangian). These three tables show that be¬ 
tween 16® X 256 and 20® x 256, there is little lattice vol¬ 
ume dependence in Mis ^nd |i?i 5 ( 0 )|®/M®g but between 
16® X 256 and 12® X 256, there are some lattice volume 
dependences. For Ma > 4.0, results from two larger lat¬ 
tice volumes agree within error bar. For Ma < 4.0, Mis 
shows no difference and |i?is(0)p/M®g. has ~ 3 % dif¬ 
ference between the two larger lattices. Thus, we can 
conclude that the finite volume effect is small in the re¬ 
sult from 16® X 256 lattices. 

As discussed, in this work stop mass is not tuned and 
a mean field estimate of the zero energy shift is used. 
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FIG. 2. The effective mass plot (bottom) and the Ao plot 
(top) for the S-wave channel with 0{v^) lagrangian using 
neighboring two points of the correlators (G(ri), G(Ti+i)) on 
16^ X 256 lattices. Error bar is from Jackknife analysis 



FIG. 3. |i?is(0)p/Mf3 from nonrelativistic correlators for the 
S-wave channel on a 16® x 256 lattice with 0{v^) lagrangian 


These two aspects are related to each other since tuning 
amounts to changing lattice stop mass until the kinetic 
mass in the nonrelativistic dispersion relation is equal to 
the mass of the a given state where part of the “correct 
mass” is from the zero energy shift, 2{ZmM — Eq). How¬ 
ever, tuning stop mass imprecisely is not a big source 
of systematic error. For example, in a bottomonium 
study |2^, the difference between a properly tuned bot¬ 
tom quark mass, Ma = 2.92 and a rough estimate. 


Ma 

k 

Eis a 

Ao a® 

Mis (GeV) 

|R(0)|VM®s 

1.0 

4 

0.1622(1) 

0.00483(1) 

84.29(1) 

4.22(1) X 10"® 

2.0 

2 

0.1693(1) 

0.02438(4) 

200.0(1) 

1.596(5) X 10"® 

2.5 

1 

0.1672(1) 

0.03630(6) 

250.1(1) 

1.216(4) X 10"® 

4.0 

1 

0.1553(1) 

0.1233(1) 

414.0(1) 

9.09(2) X 10-^ 

5.0 

1 

0.1455(1) 

0.2346(1) 

522.1(1) 

8.63(3) X 10-^ 

6.0 

1 

0.1343(1) 

0.4007(3) 

629.8(1) 

8.40(3) X 10-^ 

7.0 

1 

0.1220(1) 

0.6171(3) 

737.2(1) 

8.07(2) X 10"‘‘ 

8.0 

1 

0.1087(1) 

0.8617(4) 

844.4(1) 

7.49(2) X 10-^ 


TABLE III. Eis and Ao from lattice calculation (in lattice 
unit) on 20® x 256 lattices with NREFT Lagrangian. The 
result is based on 400 stoponium correlators and the error bar 
is from the jackknife analysis of the 1-exponential fit 


Ma 

k 

Eis a 

Ao a® 

Mis (GeV) 

|R(0)|VM®s 

1.0 

4 

0.1526(4) 

0.0048(4) 

83.81(2) 

4.27(39) X 10"® 

2.0 

2 

0.1635(2) 

0.0193(2) 

199.7(1) 

1.27(2) X 10"® 

2.5 

1 

0.1638(1) 

0.0310(2) 

249.9(1) 

1.040(9) X 10"® 

4.0 

1 

0.1549(1) 

0.1201(2) 

414.0(1) 

8.86(3) X 10"^ 

5.0 

1 

0.1455(1) 

0.2328(4) 

522.1(1) 

8.56(3) X 10"'‘ 

6.0 

1 

0.1344(1) 

0.3994(6) 

629.8(1) 

8.37(3) X lO"'^ 

7.0 

1 

0.1220(1) 

0.6155(8) 

737.2(1) 

8.04(3) X 10"^ 

8.0 

1 

0.1087(1) 

0.8592(9) 

844.4(1) 

7.47(3) X lO"'^ 


TABLE IV. Eis and Aq from lattice calculation (in lattice 
unit) on 12® x 256 lattices with NREFT Lagrangian. The 
result is based on 400 stoponium correlators and the error bar 
is from the jackknife analysis of the 1-exponential fit 


Ma = 2.90 (using Mj, = 4.65 GeV |1S] and the lattice 
spacing 0.1127 fm) is consistent compared with the accu¬ 
racy of our calculation. From Table |T] one observes that 
the mean field estimate for the mass renormalization ef¬ 
fect is small {Zm — 1 ~ 0.05) at ,0 = 8.751 and the varia¬ 
tion in Eqg due to change of Ma is ~ 0.02 for Ma < 4.0 
and ~ 0.005 for Ma > 4.0. Furthermore the difference 
between the leading order perturbative estimate for the 
tadpole factor uq and the Monte Carlo estimate is 1.4%. 
Thus, the mean field estimates themselves must be close 
to the perturbative estimates for each quantity. 

NREFT Lagrangian is an infinite series expansion in 
which upto 0{v‘^) is given in Eq. 0 and Eq. (|^. Each 
term other than the kinetic terirjj^q. @) comes with 
effective couplings and radiative corrections to these coef¬ 
ficients gives series expansions of as 1 4-0(0^) 

for Ci and are expected to be less than 10% since as{M) ~ 


in this case, tuning of the mass can absorb such correction 
(reparametrization invariance). 
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0(0.1) for M > 100 GeV. Similarly, the discretization er¬ 
ror due to the finite lattice spacing must be small since 
the magnitude of the improvement terms which correct 
for the finite lattice spacing effect of the spatial derivative 

Og ^^12 ^ the temporal deriva- 


2 4 

(a ° /£ 

12M / 


12M / M 

tive (a / Pi 

cive g. M 


Mav^ 


) are small. 


8fcM2 / M ' ' 8fc 

Estimating effects of the quenched approximation on 
the matrix element is difficult. In the bottomonium case, 
the comparison of S-wave function at the origin from the 
dynamical simulation na and that from the quenched 
simulation m using leading order NRQCD lagrangian 
found that the matrix element from the quenched approx¬ 
imation underestimates by ~ 40% since the distance scale 
associated with the bottomonium bound state j^) is 
larger than the scale at which the matrix elements sample 
the wave function (~ a = ^)- 

In the lattice calculation of the matrix elements, the 
factorization scale is related to the lattice cutoff and the 
effective cutoff is affected by the specific form of the lat¬ 
tice action and the evolution equation of lattice Green’s 
functions m- Such an effect needs to be studied if we 
are interested in the lattice matrix elements beyond the 
leading order of nonrelativistic expansion. Also, since 
the matrix element in Table |TT] is actually a lattice ma¬ 
trix element, one needs to calculate perturbative match¬ 
ing coefficients between the lattice regularization scheme 
and a continuum regularization scheme which is used for 
the parton level cross section (e.g., MS) to obtain con¬ 
tinuum matrix elements. In this leading order NREFT 
Lagrangian study, the matching is not performed. How¬ 
ever, since is small, we expect that renormalization 
effect will be small. 


IV. SUMMARY 


Under the assumption that stop is long-lived enough 
to form bound states, stoponium plays an important role 
in the study of stop and offers an interesting probe to 
stop searches mil]. In this case, nonperturbative quan¬ 
tity, \Ris{0)\'^/Mfg where Mis is the mass of l^ state of 
stoponium and i?is(0)| is the radial part of the wavefunc- 
tion at the origin, naturally appear in the study of pro¬ 
ductions and decays of stoponium. Thus far, a potential 
model estimate for \Ris{0)\'^/M^g is used in phenomeno¬ 
logical investigations of stoponium, which is unsuitable 
for improved perturbative calculations. 

In this work, using lattice nonrelativistic formulation 
for heavy stop which is interacting through the strong in¬ 
teraction, we calculate the mass of S-wave stoponium and 
the stoponium matrix element relevant for stop-anti stop 
annihilation decay, where this matrix element is related 
to |i?i5(0)p in NREFT scheme. Gompared to potential 
model studies, lattice study of stoponium is advantageous 


in that a particular functional form of the potential be¬ 
tween stop and anti stop needs not be assumed and errors 
associated with a lattice calculation can be systematically 
studied and improved. 

Monte Garlo samples of SU{3) color gauge field at 
/3 = 8.751 (a“^ = 50GeV) on 16^ x 256 lattices (which 
is generated in quenched approximation with multi-hit 
Metropolis/over-relaxation algorithm) is used for lattice 
calculation. Table |TT] and Fig. summarizes the lat¬ 
tice result for the IS” stoponium mass and the lattice 
matrix element for \Ris{0)\‘^/M^g with 0(v'^) NREFT 
Lagrangian. In general, the lattice \Ris{0)\'^/Mfg is fac¬ 
tor 3.5 ^ 4 larger than a potential model estimate of 
[TH] for 200 < Mis < SOOGeV. The difference be¬ 
tween the lattice result and the potential model result 
is larger at lighter Mis and becomes smaller at heavier 
Mis- According to this trend, the lattice result may ap¬ 
proach Goulombic behavior of the wavefunction at heav¬ 
ier Mis > ITeV although upto ~ 800 GeV the lattice 
result is still far from reaching Goulombic limit. Eis, 
the energy of l^ state of stoponium ranges from ^ 6 
GeV {Ma = 1) to ~ 8 GeV {Ma = 8). Upto order, 
toponium is equivalent to stoponium except the spin de¬ 
generacy and the result in Table is equally applicable 
to toponium since the spin-flip term in NRQGD is 0{v'^) 

m- 

Although further studies on systematic effects in 
our lattice result is necessary, factor 3.5 ^ 4 larger 
\Ris(0)\'^/Mfg implies an enhanced stoponium produc¬ 
tion cross section in hadron colliders such as LHG and can 
have interesting consequences for stop search in hadron 
colliders. Therefore, we hope to study stoponium using 
SU{3) lattice gauge fields which includes vacuum polar¬ 
ization effect of light dynamical quark and improve the 
above result using 0{v^) NREFT Lagrangian in the fu¬ 
ture. For a next-to-leading order NREFT Lagrangian 
study, other matrix elements which are higher order in 
is also needed and the perturbative matching men¬ 
tioned in the above needs to be performed. Since excited 
states of S-wave stoponium will contribute to the pro¬ 
duction cross section, matrix elements for excited states 
will be also interesting to study. 
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